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The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order 
Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au+Au 
collisions at = 200 GeV for various collision centralities. Combining two different analysis 

techniques, results were obtained in the transverse momentum range of 0.4 < pr < 4.0 GeV/c. At 
low pt the second-order coefficients, V 2 , are similar to the ones observed in hadrons. Third order 
coefficients, V3, are nonzero and almost independent of centrality. These new results on V2 and 
V3, combined with previously published results on yields, are compared to model calculations that 
provide yields and asymmetries in the same framework. Those models are challenged to explain 
simultaneously the observed large yield and large azimuthal anisotropies. 
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I. INTRODUCTION 

Direct photons emerging from relativistic heavy ion 
collisions have long been considered an important probe 
of the entire evolution of the colliding system [J. At 
almost all known or conjectured stages of the collision 
there are processes producing photons. Unlike hadronic 
observables that mostly encode the state of the medium 
at freeze-out, photons are emitted at all times throughout 
the rapid evolution of the heavy ion collision and leave 
the interaction region unmodified. Thus by measuring 
direct photons one has access to information about the 
properties and dynamics of the medium integrated over 
space and time. The measurement of direct photons is 
challenging due to a large background of photons from 
the vacuum decay of final state hadrons (tt^, tj, uj, etc.). 

The PHENIX experiment at the Relativistic Heaw 
Ion Collider reported large direct photon yields 
with strong centrality dependence Q and significant az¬ 
imuthal anisotropy or “elliptic flow” Q. Particularly 
surprising is the discovery of large azimuthal anisotropy 
for direct photons which is comparable to that ob¬ 
served for hadrons [ 3 . Preliminary results from the 
Large Hadron Collider ill indicate similar direct pho¬ 
ton yields and anisotropies. The observation of large 
azimuthal anisotropy combined with observations pub¬ 
lished earlier that the direct photon yields themselves 
are large i i contradicts several existing interpreta¬ 
tions where the large yields are provided at the very 
early production stage, when the temperature of the sys¬ 
tem is highest but the collective flow including azimuthal 
asymmetry is negligible. Conversely, the observed large 
anisotropy suggests that photon production occurs at 
very late stages of the collision when the collective flow 
of the system is fully developed, while the temperature 
and the corresponding thermal photon emission rates are 
already lower. Indeed, theoretical models have great diffi¬ 
culty to simultaneously describe the observed yields and 
anisotropy. This failure, colloquially called “the direct 
photon puzzle”, triggered a large amount of theoretical 
work, new models and insights [sl-l^. 

In this paper we present new, more precise results on 
the azimuthal anisotropy of direct photon emission from 
200 GeV Au-I-Au collisions recorded in 2007 and 2010 by 
the PHENIX experiment. Results include second and 
third order Fourier components of azimuthal distribu¬ 
tions (u 2 and U 3 , respectively) measured over a transverse 
momentum range extended down to 0.4 GeV/c. The new 
data, together with published results on yields, are com¬ 
pared to some of the more recent model calculations. 
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The paper is organized as follows. In Sec. [TTl we de¬ 
scribe the experiment, the data set, the way events are 
selected and categorized, and the two methods by which 
photons are measured. In Sec. Hill the steps needed to 
determine the direct photon V 2 , V 3 , and their uncertain¬ 
ties are described, and the final results are presented. In 
Sec.|IV]the results are compared to a few models treating 
yields and azimuthal asymmetries in a consistent frame¬ 
work. Sec. El summarizes our findings. 


II. EXPERIMENTAL SETUP AND PHOTON 
MEASUREMENTS 

In PHENIX photons are detected by two substantially 
different techniques. The first technique uses external 
conversion of photons as described in detail in Ref. 0. 
This method provides a high purity photon sample with 
good momentum resolution, but requires large statistics 
due to the few percent conversion probability and re¬ 
duced acceptance. Therefore the pr range is limited. 
The second technique is a traditional calorimetric mea¬ 
surement of photons similar to Ref. Q, but with higher 
statistics. For photons identified by either technique, the 
azimuthal anisotropy is extracted with the event plane 
(EP) method. Here we give a brief summary of the 
PHENIX detector systems and a short description of the 
two analyses. 


A. Event selection and centrality determination 

Data from 200 GeV Au-|-Au collisions were recorded 
with a minimum bias (MB) trigger based on the signal 
in the beam-beam counters (BBG) [s^, which are located 
around the beampipe at 3.1 < \ri\ < 3.9 and cover the full 
azimuth. The minimum bias trigger requires at least two 
hits in each of the two BBGs (north and south) as well as 
a reconstructed vertex from the time-of-flight difference 
between the two sides. The efficiency of the MB trigger 
is 92.3 ± 0.4(stat) ± 1.6{sys)%. 

Collision centrality is calculated as percentiles of the 
total charge distribution in the north and south BBC. 
The centrality determination is based on percentiles of 
the total charge seen in the north and south BBC and 
takes into account small shifts in 77 coverage due to vari¬ 
ations of the collision ^-vertex. 


B. Inclusive photons via external conversion 

External conversion photons are reconstructed from 
2.6 X 10® MB = 200 GeV Au-I-Au events recorded 

during the 2010 data taking period. The event vertex in 
this dataset was \z\ < 10 cm to ensure that the magnetic 
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field is sufficiently uniform. The same sample was pre¬ 
viously used in Ref. @ to determine direct photon yield 
and its centrality dependence, where details of this anal¬ 
ysis can be found. In the rest of this paper this sample 
is referred to as “conversion photons”. 

Photons convert to e“'"e“ pairs in the readout plane of 
the Hadron Blind Detector (HBD) [s^, which is located 
at ^60 cm radial distance from the collision vertex and 
corresponds to ~3% Xq, where Xq is the radiation length. 
The electron and positron from the photon conversion 
are tracked through the PHENIX central tracking detec¬ 
tors 3^. The azimuthal direction (j) and the momentum 
p are reconstructed from the drift-chamber information, 
while the polar angle of each track is determined by a 
point measurement in the innermost pad-chamber (PCI) 
and the collision vertex. High efficiency electron identifi¬ 
cation cuts are used to reduce the hadron contamination 
in the sample. Light above a minimum threshold in the 
ring-imaging Cerenkov detector [s^ and a matching clus¬ 
ter of energy E in the electromagnetic calorimeter (EM- 
Cal) [ 3 ^ such that E > 0.15 GeV and E/p > 0.5, where 
p is the momentum, are required. The EMCal com¬ 
prises two calorimeter types: 6 sectors of lead scintillator 
sampling calorimeter (PbSc) and 2 sectors of lead glass 
Cerenkov calorimeter (PbGl). The typical energy reso¬ 
lution of the PbSc is 5E/E = 8.1 %/y/E{GeV) © 2.1 %, 
and that of the PbGl is 6 E/E = 5.9 %/-yE^Gel^©0.8 %. 


The energy resolution, just as the photon identification 
efficiency, depends on centrality and its (small) effect is 
corrected for using simulated photon showers embedded 
into real events. 

All remaining tracks with pT > 0.2GeV/c, are com¬ 
bined to pairs. Gonversion photons are identified by 
analyzing the invariant mass of the pairs. The default 
tracking in PHENIX assumes that each track originates 
at the collision vertex. Thus, if the e“'"e“ pair comes from 
a conversion of a real photon in the HBD readout plane, 
the momenta will be mis-measured and a finite mass, in 
this case about niee ^12MeV/c^, is reconstructed. Con¬ 
versely, if the momenta are re-calculated assuming the 
HBD readout plane as origin, the invariant mass is close 
to zero. Through a simultaneous cut on both mass cal¬ 
culations a sample of photon conversions with a purity 
of 99% is obtained down to pr = 0.4GeV/c The 
remaining 1% of pairs are mostly from the tt*^ Dalitz de¬ 
cays. The effect on the inclusive photon is estimated 
to be smaller than 1%. 


C. Inclusive photons and tt'^s via the calorimeter 

The PHENIX EMCal is the principal detector in the 
calorimetric analysis, which is performed in a similar way 
as in Ref. [3| • The V 2 and ua are measured simultaneously 
for inclusive photons and 7r°s. A total of 4.4 x 10® MB 
Au+Au events from the 2007 data taking period are an¬ 
alyzed. The event vertex in this sample was | 2 :| < 30 cm. 

Photon candidates in the EMCal are clusters above a 


threshold energy of 0.2 GeV that pass a shower shape cut 
as well as a charged particle veto cut by the pad cham¬ 
ber PG3 immediately in front of the EMCal. However, 
photon candidates with less than 1 GeV energy are only 
used to reconstruct tt®, but are not included in the inclu¬ 
sive pho ton sample of the calorimeter. As described in 
Ref. [37j| . the remaining hadron contamination was esti¬ 
mated by comparing GEANT simulations, verified with ac¬ 
tual data. The tt® is measured via the 27 decay channel, 
with a cut on the energy asymmetry of the two photons 

\ rp _ PI 

^ ~ Ei+e -2 ^ number of re¬ 

constructed 7 r®s is taken as the integral of the two-photon 
invariant mass distribution, with the combinatorial back¬ 
ground subtracted by the mixed event method [s^. The 
signal to background ratio at 1.0 < pt < 1.5GeV/c is 
0 . 1 , rapidly improving with increasing pT- 

For the inclusive photon measurement it is important 
to restrict the measurement to a region where the resid¬ 
ual contamination from misidentified hadrons is small. 
Therefore, in the inclusive photon sample only clusters 
with E > 1 GeV are considered. On the other hand the 
inclusive (and direct) photon results presented here have 
an upper range of 4 GeV/c, which is far from the thresh¬ 
old where two decay photons from a tt® can merge in 
the calorimeter. Within this px range a purity of larger 
than 95% is achieved. The largest contamination of the 
photon sample results from antineutrons, which are not 
removed by the charge particle veto but deposit signif¬ 
icant energy through annihilation. The systematic un¬ 
certainty from particle identification (PID) of photons is 
estimated by varying both the shower shape cut (five dif¬ 
ferent settings) and, independently, applying or omitting 
the charged particle veto cut. Results from all cut varia¬ 
tions are then fully corrected. The deviation between re¬ 
sults is 3-4%, which is quoted as systematic uncertainty 
on the inclusive photon yield. 


D. Event plane determination 

PHENIX has different detector systems to establish 
the EP, which cover different pseudorapidity ( 77 ) ranges: 
the outer and inner reaction plane detector (RxNO, 
1< I 77 I <1.5, RxNI, 1.5< |? 7 | <2. 8 ), the muon pis¬ 
ton calorimeters (MPCS, -3.7< rj <-Tl, MPGN, 3.1< 
rj <3.9), and the BBG (3.1< |ry| <3.9). All these detec¬ 
tors cover the full 27r azimuth and are sufficiently sepa¬ 
rated in r] such that we do not expect auto correlations 
between the event plane determination and the photon 
production asymmetry measured. The RxNI and RXNO 
are scintillation counter systems with a 2 cm Pb con¬ 
verter that makes them sensitive to photons in addition 
to charged particles. While these photons contribute to 
the determination of the event plane, note that they are 
separated at least Ar]=0.7 from the central region, which 
is where the photon V 2 and U 3 are measured. 

The results in this paper are obtained using the event 
planes measured by the combination of the RxNI and 
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(a) RxN Inner+Outer (S+N) 

(b) MPC (S+N) 

South : -2.8<T|<-1 

South : -3.7<r|<-3.1 

North : 1<T|<2.8 

North : 3.1 <ri<3.9 
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FIG. 1. Event plane resolution as a function of centrality for 
the RxN(I+0) detector (a) used for the final results in this 
paper, and (b) for the MFC detector used to cross-check the 
results. 


to photons. It is calculated based on V 2 , fa, and yields 
measured for charged and neutral pions; contributions 
from heavier mesons are taken into account using proper 
scaling (see Sec. IIIIBI) . As a final step (iii) the direct 
photon V 2 and are calculated statistically through a 
subtraction of the results from step (i) and (ii) weighted 
by the ratio Rj, the ratio of the yields of direct photons 
to the yield of photons from hadron decays (see Eq. El). 

A. Inclusive photon V2 and vs 

The inclusive photon V 2 and V 3 are measured with re¬ 
spect to the event plane. We employ two methods to de¬ 
termine these coefficients. For each photon the azimuthal 
angular difference {(j) — 'Ffc), with k = 2,3, is calculated. 
In the first method the coefficients are determined as the 
event ensemble average for individual bins in photon px 
and centrality: 


RxNO [ 13 . Due to the large rapidity coverage this com¬ 
bination has the best resolution. The resolution Res{'^n) 
is measured with the 2-subevent method [i^. The reso¬ 
lution for RxN and MFC is shown in Figure [U The fi¬ 
nal results are cross-checked by using the other detectors 
for the event plane determination. Despite the signifi¬ 
cant difference in resolution the measured direct photon 
anisotropies are consistent, within the systematic uncer¬ 
tainties. 


Vn = (cos {n{(j) - ■^n)})/Res{-^n)- (2) 

Here Res{'i’n) is the resolution function that accounts for 
the finite event plane resolution (see Figure EJ- 

In the second method the azimuthal distribution of 
photons in a given px and centrality bin is fitted as: 


) =^o[l + 2<cos{n((/)-^n)}], (3) 


III. DIRECT PHOTON V2 AND F3 

The photon anisotropy is measured via the coefficients 
of a Fourier decomposition of the azimuthal distributions 
of photons with respect to the event plane (4 (tI | 

^ ocl-b^[^fc„cos{n((/)-^fc)}], (1) 

where 0 is the azimuthal angle of the photon, 'I'fe is the 
orientation of the event plane for a given event, and 
Vkn are the coefficients with respect to the event 
plane. In our analysis we made and explicitly tested the 
assumption that the 2nd and 3rd order event planes are 
uncorrelated, which allows us to ignore the k ^ n terms 
and to introduce the notation V 2 and for the case k = 
n, i.e. in the rest of the paper we use V 2 = V 22 and 

■^3 = '^ 33 - 

The determination of the direct photon V 2 and V 3 pro¬ 
ceeds in three steps: (i) V 2 and V 3 are determined for 
the conversion photon sample (Section IIIBj) and for the 
calorimeter photon sample fSection FlI Cl) with respect to 
the event plane fSection IlIDI) . We refer to these coeffi¬ 
cients as inclusive photon and In the second 
step (ii), the decay photon and 'Cg®® are estimated, 
i.e. the anisotropy resulting from the decays of hadrons 


Vn =v'^/Res{^n)- (4) 

This is Eq. El for the case k = n and neglecting all 
k ^ n terms. The measured values of V 2 and V 3 (^ 21 ^ 3 ) 
need to be corrected for the event plane resolution. 

In the conversion photon method the quoted Vn val¬ 
ues come from the average cosine method, while in the 
calorimeter analysis the quoted Vn values are the average 
of the results obtained with the two methods. The dif¬ 
ference between the two methods is less than 1%. The 
results for the inclusive photon V 2 and V 3 are shown in 
Figure El Both measurements agree in the region where 
they overlap. 

B. Decay photon V2 and V3 

About 80%-90% of the inclusive photons come from 
decays of neutral mesons and exhibit an anisotropy 
with respect to the event plane that results from the 
anisotropy of the parent mesons Q. To estimate this 
contribution we use measured yields and anisotropy for 
charged and neutral pions; Vn for heavier mesons is ob¬ 
tained by KEx scaling as described below. The yields 
of mesons used here are the same as are used for the 
measurement of Rj in Ref. Q. 
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FIG. 2. Inclusive photon V 2 and V 3 at midrapidity (|? 7 | < 0.35) for Au+Au collisions at ~ 200 GeV in different centrality 

bins 0%-20% (a,d), 20%-40% (b,e), and 40%-60% (c,f) with the event plane estimated with the reaction plane detector 
(1< \r]\ <2.8). The data from the external conversion method are shown as solid circles and from the calorimeter method as 
solid squares. The error bars (boxes) around the data points are statistical (systematic) uncertainties. Also shown are the 
calculated decay photon V 2 and V 3 (thick solid line) along with the statistical (dotted line) and systematic (light solid line) 
uncertainties resulting from uncertainties on the input data. An additional systematic uncertainty due to the finite event plane 
resolution is not shown (see Tabled]), because it is common to all Vn measurements. 


The V 2 and ^3 for pions are determined by combining 
data from different measurements of charged and neutral 
pion V 2 and V 3 . The 7 r° V 2 has been published in Ref. [4l| 
but the measurement has been repeated in this analysis 
to check the consistency of the results. The method to 
count the number of 7 r°s in any px bin is briefly described 
in Sec. Ill Cl To obtain V 2 (vs) for each px the number 
of reconstructed 7 r°s is extracted in six 15 (10) degree 
wide bins of the azimuthal angle A$ = $ — dfra where 

is the azimuth of the 7 r° and 'i’n is the second (third) 
order event plane (see Sec. Ill Dl) . These distributions of 
the raw 7 r° counts vs. A$ are then fitted as described in 
Sec. IIII Al to obtain V 2 and V 3 for 7 r°. Note that because 
the individual 7 r°s are not identified, the average cosine 
method [i^ is not applicable. 

These data are combined with tt^ data in the px range 
0.5 to 4 GeV/c [i^. For V 2 we also use data from 
Ref. [ 4 ^ . For the centrality class 20%-40% these data are 
compiled in Figure [3| We interpolate the data, weighted 
by their statistical and systematic uncertainties, to ob¬ 
tain an average value for pions as a function of px- 
The result of this averaging procedure, including our es¬ 
timate of the systematic uncertainties, is also shown in 
Figure Ej 

For the heavier mesons, r],uj,p,iri , the is derived 
from the Vn of the pions by scaling with the kinetic energy 



v~[KEx) = vl{KEx), (5) 

where 


KEx = mx — m = Px + rn? — to, (6 ) 

where to is the mass of the corresponding meson. 

The yields of the heavier mesons are determined from 
the 7 r° yields at pT = 5 GeV/c using the following ra¬ 
tios: p/ttO = 0.46 ± 0.060, = 0.83 ± 0.12, p/tt^ = 

1.00 ± 0.300 and p /t:° = 0.25 ± 0.075. Below px = 
2GeV/c KEx-scaling is only an extrapolation for the rj 
yields. Therefore, we also applied a blast-wave fit, and 
the difference is included in the systematic uncertainties. 
Note that the blast-wave fit results in lower rj yield at 
small Px, increasing the direct photon yield and its V 2 ,V 3 . 
The meson yields, momentum spectra and v-n are used to 
simulate mesons that are then decayed to all decay chains 
including photons. From the simulation we calculate the 
decay photon using Eq. with i?es('k„) = 1, be¬ 
cause the event plane is known in the simulation. The 
only source of systematic uncertainty on u)]®® is the un¬ 
certainty of the measured 7 r° V 2 and V 3 , and the resulting 
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inclusive to decay photon yield measured in Ref. Q. 
The procedure was introduced in Ref. 0: 
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(7) 


We reproduce R^ from Ref. Q with statistical and sys¬ 
tematic uncertainties in Figured 
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FIG. 4. The inclusive over decay photon ratio used in the 
current analysis. Present data means the results published in 

Ref. 0. 
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FIG. 3. Top panels: charged and neutral pion V 2 (a) and 
Us (c) for the 20%-40% centrality class, including previously 
published results. The averaged values used in our analysis 
are shown as a thick solid line together with the estimated 
statistical (dotted line) and systematic (light solid line) un¬ 
certainties. Bottom panels: ratio of the measured V 2 (b) and 
Us (d) values to the averaged values. 


decay photon V 2 and U 3 , derived from it. The resulting 
is compared to the inclusive photon u„ in Figure [5] 
We find that the decay photon and inclusive photon v„ 
are similar. This was already observed for V 2 in Ref. [1|, 
but is now also found for U 3 . Given that a finite direct 
photon yield has already been established ii , the sim¬ 
ilarity of and implies a large direct photon V 3 , 
as will be shown in the next section. 


C. Direct photon U2 and us 

The V 2 and U 3 for direct photons are extracted from the 
measured inclusive photon u™®, the decay photon 
discussed in the previous sections, and the ratio of the 


All systematic uncertainties on the individual contri¬ 
butions on are summarized in Table HI Uncertain¬ 
ties that are uncorrelated between data points are called 
Type A, those that are correlated are Type B and un¬ 
certainties that change all points by a common multi¬ 
plicative factor are called Type C. Uncertainties on R-^ 
are common for V 2 and U 3 and for the conversion and 
calorimeter method. For photon and pion u„ measure¬ 
ments with PHENIX the orientation of the event planes, 
i.e. 'I'n, is determined with the same detectors using the 
same algorithms. Thus the systematic uncertainty on the 
event plane determination is common for all V 2 (U 3 ) mea¬ 
surements. The uncertainties on the decay photon Vn 
are common to the conversion and calorimeter method. 
The systematic uncertainty on u™® is independent for the 
two methods and mostly reflects the different purity of 
>95% compared to >99% for the calorimeter and con¬ 
version method, respectively. 

Using Gaussian error propagation, the statistical and 
systematic uncertainties would be calculated as: 
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TABLE I. Summary of systematic uncertainties on the input to the measurement of where the is from Ref. [^, and the 
and indicate “inclusive” and “decay” photons, respectively. The values are quoted for pr < 3 GeV/c, although most 
do not vary with pr, as can be seen from Figures [2] and [3] The uncertainties on the due to the statistical uncertainty of the 
input data are uncorrelated between data points (type A); they are included in the statistical errors on the final results. Type 
B uncertainties are correlated in pr, i.e. they can vary with pr but only smoothly in the quoted range. Type C uncertainties 
change for all pr by a constant multiplicative factor. The systematic uncertainties on V 2 and V 3 are typical values. 


Input 

Source 

0%-20% 

centralities 

20%-40% 

40%-60% 

Type 

Ry 


5.5% 

5.5% 

5.5% 

B 

vT 

conversion method 

<1% 

<1% 

<1% 

B 


calorimeter method 

4% 

3% 

4% 

B 


meson V 2 (stat) 

<1% 

<1% 

<1% 

A 


7r° V 2 (sys) 

5% 

3% 

2% 

B 


T],uj V 2 (sys) 

<1% 

<1% 

<1% 

B 


Event plane 

3% 

3% 

3% 

C 

vr 

conversion method 

<1% 

<1% 

<1% 

B 


calorimeter method 

5% 

7% 

10% 

B 


meson V 3 (stat) 

1% 

2% 

4% 

A 


7r° V 3 (sys) 

11% 

11% 

11% 

B 


ri,oj V 3 (sys) 

~ 1% 

~1% 

~1% 

B 


Event plane 

6% 

7% 

18% 

C 


Except for the case there is a nonlinear 

dependence on Rj that, combined with uncertainties of 
20%-30% on {Rj — 1), results in asymmetric uncertain¬ 
ties, which are not described by Eq. [S] In particular, for 
the case the uncertainties on and are 

amplified if Rj is small. 

We estimate these asymmetric uncertainties by mod¬ 
eling a probability distribution for possible values of 
using the statistical and systematic uncertainties on 
Vn^, R^, and the event plane resolution. We assume 
that the individual statistical and systematic uncertain¬ 
ties follow Gaussian probability distributions. The prob¬ 
ability distribution for is then determined by gener¬ 
ating many combinations of u™®, and R^. Figure [5] 
shows one example of a probability distribution based on 
the systematic uncertainties on the calorimeter measure¬ 
ment for 0%-20% centrality and 1 < pp < 1.5GeV/c. In 
Figure [S] the effect of the uncertainty of only u(^®® or 
R-y, are plotted separately. The asymmetry due to the 
uncertainty of Ry is clearly visible. 

Probability distributions based on statistical (includ¬ 
ing type A systematics) and systematic uncertainties are 
determined for each data point in px and centrality 
and for both analyses. The central value for each data 
point was calculated using Eq. [T] We note that the peak 
or median of the probability distributions used to deter¬ 
mine the statistical and systematic uncertainties agrees 
with the calculated central value to better than the sym¬ 


bol size. From each distribution we calculate the lower 
and upper bound on the uncertainty by integrating from 
±00 to a Vn for which the integrated probability reaches 
15.9%. These values bracket a 68 % probability range 
for Vn and are quoted as upper and lower statistical and 
systematic uncertainties on the final result. 

The final results for the direct photon V 2 and 1 ) 3 , in¬ 
cluding statistical and systematic uncertainties as out¬ 
lined above, are shown in Figure [5] for three centralities 
and separately for the two analysis methods. For the con¬ 
version method V 3 is shown only for the highest centrality 
bin; the statistical fluctuations preclude any meaningful 
measurement in the more peripheral bins. The data and 
their uncertainties are shown in Tables m and mu 

The two analysis techniques are very different but the 
results agree well in the overlap region, and they are also 
consistent with the results published earlier [^. The di¬ 
rect photon V 2 centrality dependence, both in trend and 
magnitude, is quite similar to the observed pion V 2 - The 
third order coefficients ^3 are consistent with no central¬ 
ity dependence. 

IV. COMPARISONS TO MODELS 

As already mentioned, the essence of the “direct pho¬ 
ton puzzle” is that current theoretical scenarios have dif¬ 
ficulties explaining the large direct photon yield and az- 
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TABLE II. Direct photon V 2 for the indicated centrality bins for the two methods used. Uncertainties are shown separately as 
upper and lower. 


Centrality 

Method 

< Pt > [GeV/c ] 

V2 

Statistical uncert. 

Systematic uncert. 

0%-20% 

conversion photon 

0.50 

0.0531 

4-0.0084, -0.0076 

+0.0200, -0.0187 



0.70 

0.0387 

+0.0070, -0.0087 

+0.0252, -0.0291 



0.90 

0.0357 

+0.0080, -0.0104 

+0.0185, -0.0246 



1.10 

0.0456 

+0.0105, -0.0135 

+0.0208, -0.0277 



1.30 

0.0713 

+0.0116, -0.0128 

+0.0185, -0.0207 



1.50 

0.0979 

+0.0162, -0.0153 

+0.0227, -0.0214 



1.70 

0.0735 

+0.0148, -0.0160 

+0.0157, -0.0173 



1.90 

0.1560 

+0.0291, -0.0229 

+0.0254, -0.0192 



2.25 

0.1034 

+0.0247, -0.0243 

+0.0223, -0.0215 



3.00 

0.0699 

+0.0316, -0.0338 

+0.0140, -0.0155 



4.25 

-0.3534 

+0.8077, -0.1197 

+0.1149, -0.1831 


calorimeter 

1.19 

0.0591 

+0.0038, -0.0058 

+0.0225, -0.0266 



1.69 

0.0852 

+0.0029, -0.0035 

+0.0163, -0.0170 



2.20 

0.0957 

+0.0046, -0.0050 

+0.0214, -0.0218 



2.70 

0.0903 

+0.0074, -0.0078 

+0.0186, -0.0190 



3.20 

0.0747 

+0.0098, -0.0122 

+0.0177, -0.0189 



3.85 

0.0339 

+0.0282, -0.0430 

+0.0218, -0.0298 

20%-40% 

conversion photon 

0.50 

0.0964 

+0.0125, -0.0113 

+0.0133, -0.0113 



0.70 

0.0668 

+0.0173, -0.0289 

+0.0336, -0.0485 



0.90 

0.0640 

+0.0178, -0.0281 

+0.0308, -0.0555 



1.10 

0.0866 

+0.0155, -0.0217 

+0.0240, -0.0403 



1.30 

0.1251 

+0.0146, -0.0170 

+0.0178, -0.0240 



1.50 

0.1405 

+0.0182, -0.0202 

+0.0185, -0.0227 



1.70 

0.2074 

+0.0316, -0.0269 

+0.0291, -0.0212 



1.90 

0.1511 

+0.0314, -0.0342 

+0.0207, -0.0245 



2.25 

0.1846 

+0.0279, -0.0273 

+0.0186, -0.0174 



3.00 

0.1412 

+0.0407, -0.0431 

+0.0137, -0.0160 



4.25 

0.1561 

+0.1048, -0.0992 

+0.0133, -0.0121 


calorimeter 

1.19 

0.0902 

+0.0097, -0.0151 

+0.0236, -0.0377 



1.69 

0.1403 

+0.0066, -0.0104 

+0.0185, -0.0248 



2.20 

0.1649 

+0.0046, -0.0056 

+0.0188, -0.0202 



2.70 

0.1592 

+0.0071, -0.0083 

+0.0189, -0.0200 



3.20 

0.1327 

+0.0098, -0.0136 

+0.0190, -0.0216 



3.85 

0.0972 

+0.0155, -0.0277 

+0.0153, -0.0192 

40%-60% 

conversion photon 

0.50 

0.1173 

+0.0272, -0.0252 

+0.0117, -0.0086 



0.70 

0.0905 

+0.0214, -0.0266 

+0.0149, -0.0280 



0.90 

0.1128 

+0.0261, -0.0327 

+0.0192, -0.0349 



1.10 

0.1101 

+0.0338, -0.0444 

+0.0243, -0.0473 



1.30 

0.1978 

+0.0325, -0.0313 

+0.0163, -0.0138 



1.50 

0.1608 

+0.0465, -0.0508 

+0.0168, -0.0244 



1.70 

0.3154 

+0.0943, -0.0687 

+0.0771, -0.0366 



1.90 

0.1848 

+0.0943, -0.0969 

+0.0184, -0.0224 



2.25 

0.0173 

+0.1036, -0.1478 

+0.0584, -0.1188 



3.00 

0.2305 

+0.2262, -0.1954 

+0.0473, -0.0310 



4.25 

-0.0043 

+0.4198, -0.2826 

+0.0466, -0.0920 


calorimeter 

1.19 

0.0960 

+0.0147, -0.0247 

+0.0226, -0.0462 



1.69 

0.1412 

+0.0139, -0.0255 

+0.0162, -0.0334 



2.20 

0.1687 

+0.0172, -0.0258 

+0.0212, -0.0313 



2.70 

0.1624 

+0.0323, -0.0427 

+0.0302, -0.0405 



3.20 

0.1388 

+0.0539, -0.0657 

+0.0319, -0.0487 



3.85 

0.0999 

+0.0670, -0.0788 

+0.0290, -0.0533 

imuthal asymmetries at the same time. 

This is illustrated 

Over the past few years many new 

ideas have been pro- 


by a recent state-of-the-art calculation of viscous hydro- 
dynamic calculation of photon emission with fluctuating 
initial density profiles and standard thermal rates [13, 
which falls significantly short in describing yield and V 2 . 


posed to resolve this puzzle, including non equilibrium 
effects [3 Si EE E^. e nhanced early emission due to 
large magnetic fields [3 EE E3 , enhanced emission at 
hadronization 1311 , as well as modifications of the forma- 
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TABLE III. Direct photon V 3 for the indicated centrality bins for the two methods used. Uncertainties are shown separately 
as upper and lower. 


Centrality 

Method 

< Pt > [GeV/c ] 

^3 

Statistical uncert. 

Systematic uncert. 

0%-20% 

conversion photon 

0.50 

0.0094 

-h0.0155, -0.0163 

+0.0039, -0.0052 



0.70 

0.0237 

+0.0142, -0.0146 

+0.0099, -0.0111 



0.90 

0.0094 

+0.0143, -0.0163 

+0.0119, -0.0173 



1.10 

0.0333 

+0.0204, -0.0218 

+0.0193, -0.0223 



1.30 

0.0558 

+0.0247, -0.0247 

+0.0233, -0.0233 



1.50 

0.0299 

+0.0314, -0.0346 

+0.0246, -0.0301 



1.70 

0.0476 

+0.0305, -0.0317 

+0.0161, -0.0177 



1.90 

-0.0006 

+0.0461, -0.0535 

+0.0189, -0.0265 



2.25 

0.2094 

+0.0657, -0.0516 

+0.0461, -0.0299 



3.00 

0.0637 

+0.0672, -0.0672 

+0.0172, -0.0174 



4.25 

0.2753 

+0.4140, -0.4118 

+0.1492, -0.0765 


calorimeter 

1.19 

0.0298 

+0.0055, -0.0073 

+0.0214, -0.0256 



1.69 

0.0461 

+0.0040, -0.0053 

+0.0166, -0.0182 



2.20 

0.0587 

+0.0096, -0.0110 

+0.0170, -0.0185 



2.70 

0.0696 

+0.0129, -0.0129 

+0.0180, -0.0180 



3.20 

0.0726 

+0.0191, -0.0175 

+0.0231, -0.0221 



3.85 

0.0677 

+0.0380, -0.0332 

+0.0408, -0.0378 

20%-40% 

calorimeter 

1.19 

0.0178 

+0.0085, -0.0127 

+0.0240, -0.0343 



1.69 

0.0415 

+0.0108, -0.0154 

+0.0304, -0.0381 



2.20 

0.0619 

+0.0128, -0.0146 

+0.0339, -0.0365 



2.70 

0.0703 

+0.0198, -0.0206 

+0.0326, -0.0336 



3.20 

0.0637 

+0.0244, -0.0256 

+0.0274, -0.0284 



3.85 

0.0308 

+0.0265, -0.0331 

+0.0228, -0.0250 

40%-60% 

calorimeter 

1.19 

0.0346 

+0.0131, -0.0157 

+0.0362, -0.0422 



1.69 

0.0638 

+0.0271, -0.0273 

+0.0497, -0.0494 



2.20 

0.0920 

+0.0651, -0.0567 

+0.0780, -0.0676 



2.70 

0.1011 

+0.1224, -0.1028 

+0.0973, -0.0793 



3.20 

0.0187 

+0.1580, -0.1476 

+0.0823, -0.0877 



3.85 

-0.0430 

+0.1421, -0.1289 

+0.0751, -0.0938 


tion time and initial conditions [ 2 ^ [13, [13| . 

In this subsection we compare our results to a subset 
of the models which (i) consider thermal radiation from 
the QGP and HG (hadron gas) plus additional proposed 
sources, (ii) have a complete model for the space-time 
evolution, and (iii) calculate absolute yields and V 2 . For 
the comparison we use the data for the 20%-40% central¬ 
ity class, and note that the comparison leads to similar 
conclusions for the other centrality bins. While none of 
the models describes all aspects of the available data, 
they are representative of how different theories are try¬ 
ing to cope with the challenge. 

First, we compare the data to the “fireball” scenario 
originally calculated in Ref. [11]. The model includes 
pQCD, QGP and HG contributions, with the instanta¬ 
neous rates convoluted with a fireball expansion profile. 
The basic parameter is the initial transverse acceleration 
of the fireball, ot- The prompt photon component is esti¬ 
mated in two ways. The first variant is a parametrization 
of the photon yields measured in p+p by the PHENIX 
experiment [i^ (labeled as “primordial 1 ”), the second is 
an xt-scaling motivated parametrization (labeled as “pri¬ 
mordial 2”), modified with an empirical factor K = 2.5 
to match the measured data at high px (above 4 GeV/c). 
The yield calculation includes thermal yields from the 


QGP with To = 350 MeV and from the hadronic phase. 
Different from an earlier version of the model, chemical 
equilibrium prior to kinetic freeze-out is no longer as¬ 
sumed. This results in a large enhancement in photon 
production in the later hadronic stages via processes like 
meson annihilation (for instance tt -|- p —>■ tt -|- 7 ). With 
an initial transverse acceleration ax = 0 . 12 c^/fm and 
T Ri 15fm/c fireball lifetime, 100 MeV freeze-out temper¬ 
ature and /3s = 0.77 surface velocity, the observed low 
px photon yields are recovered within systematic uncer¬ 
tainties, but underpredict the data [ll|. In Figure [7] the 
data are compared to the most recent updated “fireball” 
scenario shown in Ref. [l^, which includes a calculation 
with ideal hydrodynamics with finite initial flow at ther- 
malization and enhanced yields around chemical freeze- 
out temperature Tc that improves the description of the 
data. The direct photon V 2 has its maximum at about 
the same px in both theory and data. The V 2 calcu¬ 
lated in the original fireball scenario [ll| under predicts 
the measured one. The radial boost hardens the photons 
from the hadronic gas (HG) and in this way increases V 2 
as well, but the calculation still falls short of the mea¬ 
surement. V 3 is currently not calculated in this model. 

Second, in Figure [5] the data are compared to three 
calculations evaluated with the hydrodynamical back- 
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Direct photon 


FIG. 5. This example shows the direct photon V 2 ^'^ measured 
via the calorimeter method with the event plane estimated 
by the reaction plane detector (1 < I 77 I < 2 . 8 ) in the 0 %- 20 % 
centrality bin. Each of the various dashed curves indicate the 
probability distribution of the result due to the variation 
of a single term in Eq. [7] While varying and alone 
leaves the uncertainty on 112 '^ Gaussian, varying R-y results in 
strongly asymmetric shapes. The black solid curve shows the 
result when all uncertainties are taken into account simulta¬ 
neously. 

ground as described in Ref. [i^, The first calcu¬ 
lation, labeled ”QGP w/ viscous”, was evaluated using 
the AMY photon emission rate in the high-temperature 
(QGP) region, and included viscous corrections to the 
photon emission rates due to both bulk and shear 

viscosities. The same calculation without the viscous cor¬ 
rections corresponds to the curve labeled ”QGP, w/o vis¬ 
cous” . Once viscous corrections are included, V 2 drops by 
more than 50% at 3 GeV/c, while the yield decreases just 
by ^10%. The third curve, labelled ”semi-QGP, w/o vis¬ 
cous”, shows the consequence of including the effect of 
confinement on the photon emission rate, as computed 
in the semi-QGP approach [l^. The utilization of the 
semi-QGP photon rates at high temperatures suppresses 
the spectrum, but does not change the significantly. 
This is a consequence of the small contribution of QGP 
photons to the thermal photon V 2 , which is dominantly 
produced at temperatures around and smaller than the 
confinement temperature. The prompt photon contribu¬ 
tions in all three calculations are evaluated within the 
perturbative QGD framework. 

Third, we compare the data with PHSD (parton- 
hadron-string dynamics), a microscopic transport 
model [T^. In addition to the traditional QGP and HG 
sources (resonance decays) this model includes late stage 
meson-meson and meson-baryon Bremsstrahlung, which 
enhances the yield at the lowest pt substantially and in¬ 
creases V 2 by almost 50% in the pt < 3 GeV jc region (see 
Figure 2 in Ref. El)- Gontributions from photonic de¬ 
cays of (p and oi are also included, because these are not 


subtracted in the measurement. After all other sources 
are added, the direct photon spectrum is very well re¬ 
produced below 3 GeV/c, but V 2 under predicts the mea¬ 
sured values. Also, the pt where V 2 reaches its maxi¬ 
mum is under predicted. In Figure [9] the data are com¬ 
pared to the latest PHSD model calculation that 
included additional photon production channels in the 
hadronic phase and improved the Bremsstrahlung calcu¬ 
lation. The model also provides U 3 . It is positive and 
consistent with the data within uncertainties. 

Explaining the large yield and strong flow simulta¬ 
neously requires significant improvements in quantifying 
the contributions from the late stage QGP and hadron- 
gas interactions. Even deeper insight on both the photon 
sources and the time profile of the system may be nec¬ 
essary to further improve the models. Euture measure¬ 
ments of more differential quantities will help to distin¬ 
guish and quantify the individual photon sources. 


V. SUMMARY AND CONCLUSIONS 

The PHENIX experiment at the Relativistic Heavy Ion 
Collider measured 2nd and 3rd order Fourier coefficients 
of the azimuthal distributions of direct photons emitted 
at midrapidity in = 200 GeV Au-|-Au collisions, for 

various collision centralities. Two different and indepen¬ 
dent analyses are used to determine the inclusive photon 
yield. The external conversion photon measurement al¬ 
lows one to extend the pr range down to 0.4 GeV/c com¬ 
pared to 1.0 GeV/c for the calorimetric measurement. In 
the overlap region the two results are consistent. The V 2 
measurements are also consistent with earlier published 
results, while U 3 is published for the first time. 

Both the direct photon V 2 and ^3 are found to be large. 
The V 2 exhibits a clear centrality-dependence, while V 3 is 
consistent with no centrality dependence. At all central¬ 
ities, the direct photon V 2 is similar in magnitude to the 
hadron V 2 for px < 3 GeV/c, The direct photon V 3 is con¬ 
sistent with that for hadrons over the entire px range. 

We compare the data to several recent calculations, 
which treat the direct photon yields and the azimuthal 
asymmetries in a consistent production and evolution 
framework. None of them describe the full systematics of 
the data adequately, but there has been progress in the 
last few years. The general trend of the models appears 
to be including sources from the earliest (pre-equilibrium, 
see for instance Ref. El) or very late times in the evo¬ 
lution of the system, while giving less emphasis to pho¬ 
ton production at intermediate times, when most of the 
expansion occurs. PHSD includes new sources from the 
hadron gas and photon production even after the hadrons 
are decoupled from each other, which improves descrip¬ 
tion of the yields but still under predicts V 2 - The model 
that best approximates the measured V 2 , including the 
px region where V 2 reaches its maximum value, starts the 
evolution with a large initial boost even before thermal- 
ization El- It is also worth noting that the microscopic 
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FIG. 6. Direct photon V 2 and vs at midrapidity (|r;| < 0.35), for different centralities, measured with the conversion method 
(solid circles, green) and calorimeter method (solid squares, black). The event plane was determined with the reaction plane 
detector (1 < |7)| < 2.8). The error bars (boxes) around the data points are statistical (systematic) uncertainties. 
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FIG. 7. Comparison of the direct photon yields [3] and V 2 with 
the fireball model [I^]. The two curves for V 2 correspond to 
two different parametrizations of the prompt photon compo¬ 
nent. See text for details. 


FIG. 8. Comparison of the directphoton yields and V 2 , vz 
with a hydrodynamical model [d^, |4^ calculated under three 
different assumptions including the “semi-QGP” scenario 0. 


transport model [l3| is able to describe the anisotropies 
as well as the full-scale viscous hydrodynamics (l^ . 

While the data are getting more differential and more 
accurate, and model calculations improve, the “direct 


photon puzzle” remains unresolved. High quality data 
of yields and V 2 and for different collision systems, in¬ 
cluding very asymmetric ones, and energies would help to 
further improve our understanding of direct photon pro¬ 
duction because robust models must be able to describe 
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FIG. 9. Comparison of the direct photon yields and V 2 with 
the PHSD model 0,143]. 

the data over a wide range of experimental conditions. 
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